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ABSTRACT 

We quantify the rapid variations in X-ray brightness ("flares") from the ex- 
tremely massive colliding wind binary 77 Carinae seen during the past three or- 
bital cycles by RXTE. The observed flares tend to be shorter in duration and 
more frequent as periastron is approached, although the largest ones tend to be 
roughly constant in strength at all phases. Plausible scenarios include (1) the 
largest of multi-scale stochastic wind clumps from the LBV component entering 
and compressing the hard X-ray emitting wind-wind collision (WWC) zone, (2) 
large-scale corotating interacting regions in the LBV wind sweeping across the 
WWC zone, or (3) instabilities intrinsic to the WWC zone. The first one appears 
to be most consistent with the observations, requiring homologously expanding 
clumps as they propagate outward in the LBV wind and a turbulence-like power- 
law distribution of clumps, decreasing in number towards larger sizes, as seen in 
Wolf-Rayet winds. 

Subject headings: X-rays: stars - stars: early-type - stars: individual (77 Car) - 
stars: LBV 
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Introduction 



The extremely massive binary star rj Carinae consists of a Luminous Blue Variable 
(LBV) primary (star A) coupled with a hot, fast-wind secondary (star B) which is probably 



an evolved O star (Verner et al. 2005) or a Wolf-Rayet star (Pittard & Corcoran 2002). 



Orbital and stellar properties gleaned from the literature are given in Table [TJ These values 
are considered to be the best current values, although overall relatively uncertain. Star A 



has gone through at least one giant LBV eruption, which occurred in the 1840s (Davidson 



& Humphreys 1997). 



r\ Carinae generates X-ray emission in the 2 — 10 keV range due to the collision of star A's 
dense, slow wind with the thinner, faster wind of star B ( Ishibashi et al.|19 99). On timescales 
of years, the X-ray emission varies in a manner similar to that expected from a long-period, 
highly eccentric, colliding-wind binary, with Lx ocl/D (where D = orbital separation), for 



adiabatic conditions as is most likely the case in such a long-period system (Usov 1992) 



though near periastron passage radiative effects ( Parkin et al. [2009 ) may become important. 
This trend is broken by a broad atmospheric eclipse in the X-ray light-curve, which begins 
near periastron passage when the dense inner wind of Star A starts to block out most of the 
X-ray emission arising in the bow-shock head region of the wind-wind collision (WWC) zone. 
Along with this long-term variation, short term variations in rj Carinae's X-ray brightness, or 
"flares" ( Corcoran et al.|1997 Davidson et al.|1998 Ishibashi et al.|19~99 ), have been observed 
by RXTE for more than two full orbital cycles from 1996.2 up to at least 2009.0, near the 
date of this writing. These flares are not expected in standard colliding wind models, and 
their origin is currently unknown. We examine here the characteristics of these flares and 
discuss possible physical mechanisms for their production. 



The RXTE Observations 



The Rossi X-ray Timing Explorer (RXTE, \Bmdt et al.|1993| ) is a satellite X-ray observa- 



tory designed to provide high time-resolution, low spectral resolution observations of bright 
X-ray sources. For our purposes we restrict ourselves to discussion of data from RXTE's 
Proportional Counter Array (PCA), which consists of five collimated Proportional Counter 
Units (PCUs 0-4) sensitive to emission in the 2 — 60 keV band with a field of view of ~ 1° 
FWHM. RXTE was launched in December 1995 and began observations of rj Carinae in 
February 1996. 



The RXTE observations in the first two cycles of rj Carinae were described in Corcoran 
(2005), with similar observing parameters for the most recent data. Each PCA observation 



-3- 



Table 1. Adopted Parameters for the r] Carinae Binary System 



Parameter 



Value 



Reference 



P (orbit) 

Ta 
o 

i 
e 
to 
a 

D(periastron) 

D(apastron) 

M A 

M B 

M A 

M B 

Voo.A 
V^ B 

Ra 



2024.0 ± 2 d = 5.541 ± 0.006 yr 
HJD 2450800 ± 3 d (=1997.9618 UT) 
45° 
0.9 
243° 
15.4 AU 
1.54 AU 
29.3 AU 
90 M 
30 M 
10~ 3 M yr" 1 
10~ 5 M yr" 1 
500 km/s 
3000 km/s 
0.28 AU 



Corcoran 



Corcoran 



(2005) 



(2005) 



Okazaki et al. 



Okazaki et al. 



Okazaki et al. 



Corcoran et a' 



(2008) 



(2008) 



(2008) 



(2001) 



Hillier et al. 



(2001 



Verner et al. 



Hillier et al. 



2005) 



(2001 



Pittard & Corcoran 



Hillier et al. 



(2001) 



Pittard & Corcoran 



Hillier et al. 



(2001) 



(2002) 



(2002) 



a T D is the Heliocentric Julian Day Number at the start of the X-ray minimum in 1997-1998, 
which is likely close to the time of periastron passage. 
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lasts typically 1-2 ksec and is usually obtained at a cadence of a few times per month. For 
certain key times in the cycle (for example, just before the 2003 X-ray minimum and during 
the decline to minimum in 2008-2009) daily X-ray observations were obtained. Prior to the 
1997-1998 minimum, observations were done weekly, with daily observations only after the 
beginning of the X-ray minimum. Obtaining frequent observations just prior to X-ray minima 
is important to temporally resolve the flaring behavior of the source, especially as the X-ray 
maximum is approached when the flare durations tend to be much shorter. The data were 
extracted from PCU2 since it is the only Proportional Counter Unit which has continuously 
observed rj Carinae from the start of the campaign and thus the PCU2 data provide the 



most uniform set of observations available. We reduced the data as described by Corcoran 



(2005) and references therein. In addition to those steps, we also reprocessed all the data 
since 1996 to correct for newly-discovered problems in the SAA (South Atlantic Anomaly) 
history file that result in errors in background rate estimation, along with a bug in the RXTE 
background estimator (PCABACKEST) for faint models in Epoch 5c (corresponding to data 
obtained from 01 Jan 2004 to the present). 



3. The Residual Flare-Lightcurve 

Characterizing the flares is not easy, especially for most of the cycle away from the X-ray 
minimum when the flares are broad and weak. Flaring occurs on top of a slowly-varying 
phase-dependent baseline brightness-level produced by the shocked gas in the WWC, which 
(as noted above) mostly varies as 1/D, where D is the instantaneous separation between the 
centers of the two stars. In addition there can be small-scale variations in the slow baseline 
from one observation to the next because of uncertainties in the correction for instrumental 
background or, possibly, intrinsic variations in the source (or variability from another source 
in the PCA field of view). While Poisson-based errors are typically ~0.1 ct/s, background 
uncertainties can reach as high as ~1 ct/s rms. 

We define a "residual flare-lightcurve" as the height of the observed PCU2 net count rate 
above the variable-baseline colliding-wind emission in the following way. We first selected 
by eye a set of 55 data points to represent the apparent non-flare baseline brightness level 
due to the underlying colliding-wind emission. Using these data points, we then created a 
piecewise continuous function to represent the underlying (non-flare) brightness of the source 
by linearly interpolating the discrete data points to the times of the full RXTE lightcurve. 
We then subtracted this piecewise continuous curve from the background-corrected RXTE 
lightcurve to generate the "residual flare-lightcurve" . Figure [T] shows the discrete data points 
we selected along with the piecewise continuous brightness variation, and the residual "flare 
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lightcurve" after subtracting the interpolated baseline. We divide the entire residual "flare 
lightcurve" into 3 cycles by phase. Cycle 1 refers to phases —0.35 < <fi < 0.00, cycle 2 refers 
to 0.00 < < 1.00, and cycle 3 refers to 1.00 < < 2.00, where phases <fi are calculated as 

<f>=(t-T )/P, (1) 

where t is the Heliocentric Julian day number of the start of the observation and T Q and P 
are given in Table 1. Figure [2] compares the flares in a phase plot (see below for discussion). 

Flares which occur at higher X-ray baseline flux just before plummeting to X-ray min- 
imum are relatively narrow and intense and thus easy to identify. At lower baseline X-ray 
fluxes, flares are weak and broad, and determining flare parameters for them is more sub- 
jective. After selection by eye of all the most obvious flares whose peak value lies at least 
~1 ct/s above the zero level of the residual flare-lightcurve, we then measured the flare 
properties: time corresponding to the flare centroid, flare temporal full width at half maxi- 
mum (FWHM), and peak flare height above the baseline (H). Table [2] (appendix) lists the 
measured quantities for each identified flare. Figure [3] shows FWHM as a function of phase 
for each of the three cycles. Time separation for a given flare (also given in Table |2| was 
obtained by averaging the time interval from the prior flare to the current flare and the 
time interval from the current flare to the next flare. Figure [4] shows FWHM versus time 
separation. Both quantities are clearly correlated, although not exactly linearly. Figure [5] 
shows a phase plot of the total strength S of the flares, estimated as S = FWHM xH. We 
see that, within the scatter, the maximum total strength is ^independent of phase, whereas 
one sees a larger number of weaker flares closer to periastron. Clearly, all three cycles show 
the same trends with phase in all these quantities (FWHM, separation and strength). 



3.1. Comparison of Identified Flares 



Figure [2] compares the whole residual flare lightcurves in the three different cycles, 
including a close-up just before the X-ray minimum. There is no apparent detailed repetition 
of the flares from one cycle to the next, although a few of the flare phasings seem to repeat 
in two (or in rare cases, all three) cycles. Based on a limited sample of the Cycle 1 data, 
Corcoran et al. (1997) and Ishibashi et al. (1999) suggested a periodicity of ~ 85 days for the 



X-ray flares which they suggested might arise in stellar pulsations of r] Carinae A. We here 
attempt a more objective test for repeatability of the flares from one cycle to the next now 
using data from all 3 cycles observed by RXTE. For this purpose, Figure [6] first compares 
the complete residual lightcurves of cycles 1 and 3 to that of cycle 2. Because the sampling 
rates are different, for this comparison we interpolated the cycle 3 residuals to the phasing 
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of the cycle 2 data; because cycle 1 is somewhat undersampled compared to cycle 2 (and 
cycle 3), we interpolated the phasing of the cycle 2 data to the phasing of the cycle 1 data, 
then overplotted the cycle 1 residuals vs. the interpolated cycle 2 residuals. In Figure [6j 
both axes are biased to positive values (making a simple correlation analysis inappropriate), 
except for data points close to the origin, which are dominated by noise (and systematic 
uncertainty in our determination of the baseline level). 

Figure [7] shows a histogram of the residual intensities for all 3 cycles. Here, we see a 
trend of constant slope only above a critical value I c ~ 1.3 cts/s, to which we fit a straight 
line out to I = 10 cts/s (above which the bins become too sparsely populated). Below I c the 
values flatten out; we take this as the residual count-rate limit below which the data are noise 
dominated. We note that the slope of the straight-line fit in Figure [7] is b = —2.9 ±0.3. This 
slope is somewhat steeper than that found for the variable components of optical emission 



lines in WR stars: b m —2.0 ± 0.2 on average Moffat et al. (1994), where we may be seeing 



the power-law spectrum due to full-scale, supersonic compressible turbulence in their winds. 



Figure [8] shows the variation with polar angle 6 = arctan(x/?/) of the points in Fig. [6] 
which lie beyond the circle I c = (x 2 + y 2 ) 1 ^ 2 = 2 cts/s. The variables x and y refer to the 
abscissa and ordinate, respectively, of Fig. 6. The points have been binned into 10 angular 
bins summed over all intensities. Figure [8] shows no obvious dependence on angle between 
0° and 90° (which is the physically meaningful range) as would otherwise be expected for 
a positive correlation of the timing of flares in different orbital cycles. Assuming Poisson 
errors for each bin, a x 2 test of the null slope hypothesis between 0° and 90° confirms this 
impression: xl = 1-01 f° r a horizontal, straight-line fit. Thus we conclude that there is no 
statistical correlation between the phasings of flares in one cycle compared to another cycle: 
flare production is apparently a repetitive but stochastic process. 

To further search for possible periodicities in the observed flaring events in all three 
cycles, we Fourier analyzed the whole data string of residual flux using a periodogram ap- 



propriate for unevenly sampled data (Scargle 1982 Home & Baliunas 1986). The result is 



shown in Fig. [9J where the only narrow Fourier amplitude peaks (i.e. strict periodicity) 
we see are related to rj Carinae's orbital period and its numerous harmonics. Broader flat 
peaks are seen centered at periods of 67d, 25d, 15d and lid, suggestive more of the preferred 
flare time scales rather than any strict periodicity. However, the question remains whether 
the flares do in fact appear in regularly spaced angular intervals. We look at this later in 
subsection 4.1.2 in the context of one of our suggested models. 
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3.2. X-ray Hardness Variations 



The hardness ratio analysis presented in Figure 10 of Corcoran (2005) shows that some 



strong flares are associated with significant change in X-ray hardness, in the sense that 
the flare hardness increases towards the flare peak and then declines. This impression is 
confirmed by comparing an XMM X-ray spectrum from 2003 June 13 (near the peak of flare 
#23) and another spectrum from an observation on 2003 June 8, prior to the flare rise (see 



Figure 4 in Hamaguchi et al. 2007). These two spectra show that the increase in emission 



at the flare peak is confined to energies E > 3 keV for this strong flare. This suggests 
an increase in the emission measure of the hottest plasma with little change in low-energy 



absorption. Figure [10] shows a comparison of CHANDRA spectra obtained on 2008 Dec 8, 
near the peak of flare #43, along with a spectrum from an observation on 2008 Dec 12, 
during the decline of the flare. Once again most of the variation occurs at energies E > 3 
keV, with little change at lower energies. 

To crudely quantify the spectral change, we fit each spectrum with an absorbed single 



temperature thermal model using a standard interstellar absorption model (Morrison & 



McCammon 1983) plus a thermal emission plasma model (APEC, |Smith et al.1|200~Ij ) . The 
best fits to the spectra are shown in Figure 10 The best fit values for kT are kT = 
4.9±0.5 keV and kT = 3.6±0.5 keV for the peak and off-peak spectra, respectively. The best 
fit columns for the peak and off-peak spectra, respectively, are Nh = 5.0±0.2 cm -2 and Nh = 
4.9 ± 0.2 cm" 2 . Thus the flare peak spectrum can be characterized by a higher temperature 
than the off-peak spectrum, while there is no significant difference in the absorbing column 
for either spectrum. 

Whatever the mechanism which produces the flares (see §4 below), any plausible mecha- 
nism leads to a compression and thus additional heating of the shocked gas. This will in turn 
lead to an increase in temperature and thus X-ray hardness. We describe the compression 
process more explicitly in subsection 4.1.1. 



3.3. Summary of the Statistical Flare Properties 

The following is a summary of the salient features of the X-ray flares in t] Carinae: 



There is a dramatic decrease of flare time-width (FWHM) towards periastron. 

There is a dramatic increase of flare frequency (or decrease of inter-flare interval) 
towards periastron. 
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• The above two effects (time width and time interval) appear to be directly, although 
not necessarily linearly, correlated. 

• There is a dramatic increase of maximum flare peak-intensity (H) towards periastron. 

• Flares appear randomly but on preferred time scales. 

• During the strong-flare interval just before X-ray minimum, the flare spectra are harder 
than the non-flare colliding wind emission (and apparently harder than the spectra of 
flares near the middle of the cycle at apastron). 

• Within the noise, the maximum value reached for the total strength of the flares (S) is 
relatively constant with phase. However, at phases when the flares are narrower, there 
is a larger spread in observed flare strength. 

Any viable model to explain the X-ray flares in rj Carinae must take all of these factors 
into account. We now look at how to do this in a quantitative way. 



Plausible Flare Models 



The temporal behavior and spectral hardness of the flares suggest a disturbance of the 
interacting region of the two winds where X-rays arise, on random but preferred time scales. 
The flares must represent a change in X-ray emission level in the relatively small region where 
the two winds collide nearly head-on, of half opening angle Ox- This region appears to be 
significantly smaller in angular size than the overall shock cone opening angle in t] Carinae 



(cf. Henley et al. 2008). The half opening angle of the complete shock cone, 6, can be 
calculated using the ratio of wind momenta from Table 1: rj — (rri gVoo ^ / (rriAVoo, a) = 0.06 



for i] Carinae. From Figure 4 in Canto et al. (1996), this yields the asymptotic opening 
angle 9 = 42°. For such a wide system, it is unlikely that this angle will be significantly 
widened by radiative braking from the hotter star B ( Gayley et al.|[l997 ), until perhaps near 
periastron when radiative inhibition by star A of the wind from star B could be important 



(Parkin et al. 2009). In this section, we explore several plausible ways how the flares might 



occur. 



4.1. Flares from Structured Winds 



Two distinct kinds of structures (e.g. Moffat et al. 1994) are known to occur in hot- 
star winds: (a) stochastic multi-scale clumping in all winds and (b) large-scale corotating 
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interaction regions (CIRs, leading to discrete absorption components, DACs, in P Cygni wind 
line profiles) in many winds. We consider below whether the interaction of such structures 
with the X-ray emitting bow shock could conceivably produce the observed X-ray flares in 
7] Carinae. 



4-1.1. Flares from Clumped Winds 



Winds from hot stars are believed to be composed of a hierarchy of turbulent clumps, 
with few large ones and increasingly many more small ones, and with clump size and mass 



distributed as a power-law (e.g., Moffat 2008 ). Whether the driver of clumping is some kind of 



radiative wind instability (eg., Lucy 1982 Owocki et al. 1988) or due to subsurface convection 



(Cantiello et al. 2009) remains to be clarified. In particular, LBVs display the greatest 



amount of stochastic polarimetric variability among all hot luminous stars, due to the large 
number of clumps in their winds ( Davies et al.||2005 2007). One conjecture is that the X-ray 
flares in 77 Carinae are produced when the rarer, larger clumps enter into the colliding-wind 
zone, while the background X-ray emission comes from an unresolved superposition of many 
smaller clumps, which produces similar X-ray emission to that produced by the collision of 



smooth winds (Pittard 2007). This is analogous to the few discrete emission subpeaks seen 
on WR or O-star optical emission lines, versus the underlying ~constant wind component 



from thousands of small wind clumps (Lepine & Moffat 1999 Eversberg et al. 1998; Pittard 



2007 Lepine fc Moffat|2008 ). In principle the clumps could occur in either the wind from star 
A or star B, since ramming and being rammed lead to essentially the same effect. Assuming 
the clumps that produce the X-ray flares arise in the wind of star A, then the arrival of a 
big clump from star A into the shock-cone zone will push the bow shock towards star B, 
increasing the density of the shocked portion of star B's wind and thus increasing the overall 
X-ray emission and hardness for a brief time. Alternatively, the collision of a clump in the 
wind of either star with the colliding-wind bow shock might produce excess X-ray emission 
if the clump mixes into the shocked region (a very likely process for adiabatic collisions: 
Pittard 2007), though the speed of star A's wind (500 km s _1 ) corresponds to temperatures 
of only kT < 0.5 keV, much lower than the observed flare temperatures {kT ^3 — 5 keV). 

In any case, clumps should expand homologously with the wind. Near apastron the 
colliding wind shock is farther from either star than it is near periastron, and if the clump 
expands homologously with the wind this could explain qualitatively why the flares are 
narrower in time (i.e. FWHM) near periastron, i.e. clumps are smaller (and denser) at 
periastron, producing narrow sharp flares near periastron (and, conversely, weaker, longer 
flares near apastron). 
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Below we construct a simple phenomenological model of the interaction of clumps with 
the X-ray emitting portion of the bow shock. In what follows we assume that the angular size 
of the X-ray emitting portion of the bow-shock head as seen from either star is ~ constant 
for all phases. This is expected if both winds have constant mass-loss rate and constant 
expansion speed at the point of impact. Near periastron, however, the winds might collide 
at sub-terminal speeds, depending on the real periastron separation, the wind acceleration 
(usually parameterized as V(r) = V^,(l — R/r)@, where is the terminal velocity, R the 
stellar radius, r the radial distance and /3 ~ 1 for OB star winds), and on whether radiative 
effects are important. We note that recent modeling of the X-ray flux and hardness variation 



(Parkin & Pittard 2008) suggests that radiative effects may be important near periastron. 



Now we examine the radial clump sizes needed to explain the observed flare width- 



times. Figure 11 shows a sketch of this clump model. In simple colliding-wind models, 
the distance from star A to the stagnation point (where most of the X-ray emission is 
produced) is d,A = D/{\ + 77 1 / 2 ) = 0.80D while the distance from the companion (star 
B) to the stagnation point is ds = r] l ^ 2 D/{l + 77 1 / 2 ) = 0.20D. The orbital separation 
D = a(l — e 2 )/(l + ecosv), where v is the orbital true anomaly and a the semi-major axis. 
At periastron, D = a(l — e) = 1.54 AU, and = 1.23 AU, ds = 0.31 AU. At apastron, 
D = a(l+e) = 29.3 AU (19 times larger than at periastron), and = 23.4 AU and ds = 5.9 
AU. 

The shortest observed duration of a flare occurs near periastron, with FWHM « 3d 
near = 0.99 when the separation between the stars is ~ 3 AU. Assuming this to be similar 
to the time it takes for a big clump moving radially outward to penetrate the bow shock 
head and temporarily produce an increase in X-ray flux leads to a clump (radial) size of 
Ar c i < AtVoo ~ 0.9AU for star A and < 5.2AU for star B. Thus star B cannot be the source 
of the clumps since the implied clump size for star B is larger than the distance from star B 
to the stagnation point at <fi = 0.99 (~ 0.6 AU), unless the clump speed is much less than the 
terminal velocity of the wind of star B (which seems unlikely). For star A, there is sufficient 
(although not ample) room at periastron (and even more room in an absolute sense at other 
times); with a lower wind speed close to the LBV with a radius of ~ 0.3 AU, there is even 
more room. However one expects a clump moving at ~ 500 km s _1 to produce shocked gas 
with a temperature of < 10 6 K, not the k 3 x 10 7 K emission seen by RXTE during the 
flares. Nevertheless, we assume that star A is the source of the clumps, and that the excess 
density of the clump becomes mixed into the hotter shocked gas in the unstable region near 
the wind- wind collision ( Pittard||2007 ), or that as the clumps impinge on the bow shock they 



move the bow shock into denser regions of star B's hotter wind. 



What kind of X-ray luminosity enhancement can be expected when a clump traveling 
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at the terminal wind speed from star A pushes into the WWC zone? Assume that a large 
clump from star A's wind enters and dominates the WWC zone, but with density p' A > p A , 
where pa is the background density of star A's wind. Then, assuming mass conservation, 
the corresponding wind-momentum ratio within the WWC zone where X-rays arise will be 
modified temporarily to r( = r]/p, where p = p' A / Pa- This will then push the stagnation 
point closer to star B, since d B = L>r/ 1/2 /(l + v' 1/2 ) = D (v/ v) 1/2 /i 1 + (v/ v) 1/2 )- Then, 
assuming L x oc l/d B , one has L' x /L x = (p 1 ^ 2 + V 1 + f] 1 ^ 2 ) = 1-8 for the largest flares 
(as observed), assuming p = 4 as expected for strong shocks. 

Neglecting tangential motions of the bow shock and the clumps for simplicity, the flare 
time-width can be approximated by: At(FWHM) ~ Ar c i(dA)/V(dA), where Ar c i(dA) is 
the radial extent of the clump at the distance of the shock-cone head dA from star A. Then, 
homologous clump growth in the wind yields Ar c i(dA) = At^Ra^a/ Ra, where At c i(Ra) 
is the clump size near the sonic point of star A. Thus, 

At(FWHM) = k x D/(l - R'a/DY, (2) 

with constant k = Ar d (R A ) / '(V^aR'a) and r 'a = i l + V^) r a- If /3 = 0, then V = V^ A 
and the flare durations (i.e. time- widths) are (1 + e)/(l — e) = 19 x shorter at periastron 
than at apastron. In Fig. [3] we show a least squares fit of equation [2] to the observed flare 
FWHMs (assuming (3 = 1). Given the more uniform dispersion of the At(FW H M) values 
in log than in linear units, we optimized the fit in log space. 

Using the parameters from Table 1 along with rj = 0.06 and (3 = 1.0 (Hillier et al. 
2001), the best-fit clump model provides a good description of the observed FWHMs with 
k = 1.7 ± 0.1 d/AU. Using this value of k, Ar d {R A ) w 0.17 AU « 0MR A . This implies 
that the largest clumps occupy more than half the stellar radius in the radial direction near 
the sonic point, if indeed they actually start there, suggesting that the clumps are radially 
extended. At apastron the distance to the bowshock head is d A = D (apastron) /(l + y/rj) = 
23.5AU, and the clump size in the radial direction is still the same fraction of the distance 
according to homologous expansion, i.e. Ar d w 0.61d^ w 14AU, nearly equal to the adopted 
size of the semi-major axis. We also found from eq. 2 a spread in At(FWHM) from 3.4d 
at periastron to 50d at apastron. Since the time of periastron passage T Q was assumed to 
coincide with the beginning of the X-ray eclipse, we refitted eq. [2] to the data allowing T Q to 
vary. We found that a slightly earlier value of T by 7 ± 4d fits but this fit is only marginally 
better. Since there is little evidence at other wavelengths that periastron occurs before the 
X-ray minimum, we therefore ignore this and retain the assumption that periastron passage 
occurs at the start of the X-ray minimum. 

On the other hand, assuming that (1) clumps form stochastically but on a quasi- 
continuous and constant basis per unit time and in all directions around rj Carinae, (2) 
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they expand truly homologously in all three dimensions (radial and two transverse), and (3) 
they collide in a constant colliding-wind zone angle, there is no reason to expect the sepa- 
ration time between successive flares, At(sep), to vary around the orbit. This is contrary 
to what is seen in Fig. |4j however, where the variation in average flare time-separation is 
clearly correlated with the variation in FWHM. 

The correlation between the observed flare FWHMs and true separation may be a 
consequence of reduced detectability of clumps towards apastron. Quantitatively, this can 
be explained as follows. Flares are only detected when the residual flare height H > H c « 
1 ct s _1 , the threshold value derived above. We further assume that the total strength of 
a flare S follows a power law, N(S)dS = N S~ a dS, where a and N a are constants. This 
is inspired by the observed power law in Fig. 7 and similar behavior for sub-peak spectral 
intensities from clumps in WR winds 
W = FWHM takes on a preferred value, we can only detect flares with S > WH C . This 
leads to the total number of flares detected per unit time 

pSmax AT 

N tot (W) = / N(S)dS = j-^-AS l ^ x - (WH c ) l -% 
Jwh c I 1 - <*) 

where S max is the maximum value of S from the largest clumps. (Note that once absorbed 
into the adiabatic WWC zone, hydrodynamic models predict that the clumps will be quickly 
destroyed, i.e. merged into the ambient material: Pittard 2007.) Then the mean separation 
between successive detected flares is 

At(sep) = l/N tot (W) = hW^/ih - W*- 1 ), 

with constants k\ = (S max / 'if c ) Q-1 and k 2 = (a — 1)S^^/N . Reversing this relation, we 
can also write 

W = [k 1 At(sep)/(k 2 + At(sep))] 1/(Q - 1) . (3) 

A log-log fit of this last relation (with W and At in days; more mathematically stable 
than the above equivalent) is shown in Fig. 4, with \ 2 ~ minimized parameters a = 1.7 ± 
0.2, k\ = 23 ±24, k 2 = 69 ±52. While k\ and k 2 are poorly constrained (although necessary to 
characterize the relation), the power-law index a is very close to that found for WR clumps: 
a = 2.0 ± 0.2 (Moffat 1994). From Figure |H] we find that S max w 10 7 counts, and, with 
^i/(a-i) _ yields n c ~ 10 7 /A;^ <Q ~ 1 ' ) = 1.3 cts/sec, very close to that above. 

Overall, the fact that the clump model fits so well in both At(FWHM) and in At(sep) 
implies that the actual rate of clump production must indeed be uniform with time and 
the clumps do survive at least from 1.2 to 23 AU, the separation of the bow-shock head 
from periastron to apastron. We can now calculate the total number of observable flares 



Moffatfl994 ). Then at a given orbital phase 0, where 
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at periastron and extrapolate to the whole wind to estimate the total number of clumps. 
Taking W = 3d near periastron, we find N tot (W) = {k x - W°- 7 )/{k 2 W - 7 ) « 0.14 flares/d, 
or 0.49 flares/3. 5d, i.e. during one wind-expansion time Ra/Voq. If one flare corresponds to 
one clump of angular diameter A6° and clumps are distributed uniformly in a spherically 
symmetric wind around the star, then we find 0.49 x 4.1 x 10 4 /(A#°) 2 = 2300, 205, 23 clumps 
per wind expansion-time for A8° = 3, 10 and 30, respectively. If A9 = 3°, as in WR winds 
(Dessart & Owocki 2002 but see Dessart fc Owocki||2003 ) then the first value is most likely, 
while Davies et al. (2007) find A8° = 12° for LBV winds. These are compatible with the 
total number (~ 10 3 ) of clumps found in the inner parts of LBV winds by Davies et al. 
(2007) based on polarimetry. 



4-1.2. Flares from Large-scale Co- Rotating Interaction Regions 

Winds from hot stars are also widely observed to have rather stable structures (spiral 
density waves or "co-rotating interaction regions", CIRs, Mullan||1986 ) which are apparently 



tied in some way to the stellar surface and which are swept around with the star as it 
rotates. Such structures produce "Discrete Absorption Components" (DACs) observable in 
the absorption troughs of unsaturated P Cygni wind lines. One possibility for the production 
of the X-ray flares is that the orbital motion of the companion carries it through such a stable 
density enhancement in the wind, which is itself swept around by the rotation of the star 
(see sketch in Fig. 12). The mechanical effect of this interaction would be the same as for 
clumps: this interaction of the bow shock with a stable density enhancement could push the 
bow-shock towards the companion, to regions of higher density in the companion's wind, thus 
increasing the emission measure of the hot shocked gas. The emission measure enhancement 
would only last as long as it takes for the relative motion of the companion to carry it 
through the CIR. 

If the detectability of CIR passages suffers from the same problem as that for the 
clumps (i.e. slower sweeping at apastron, making for lower-intensity, broader flares, many of 
which are more likely to remain below the noise threshold), then the rate of sweep must be 
determined when they are best seen, i.e. at periastron, where At(sep) <10 d (cf. Table 2; 
if detectability is still a problem even at periastron, then this limit could be decreased even 
more). Then for n CIRs distributed more or less evenly around star A's equatorial plane 
(which is assumed to be aligned with the orbital plane for simplicity), the rate of sweep will 



be d6cm/dt > (27r/n)/10 rad/d. For the most plausible number of CIRs, n = 2 (Kaper 



et al. 1999) this yields d6 C m/dt > 114 rad/yr. 



Now, the net sweep rate of CIRs as seen by the companion (and the shock zone) is 
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d9cm/dt = dv/dt — Q, where dv/dt is the sweep rate due to orbital motion (varying from 
0.137 rad/yr at apastron to 49.4 rad/yr at periastron) given by 

dv 2n (1 + e cosv) 2 
dt = ~P (1-e 2 ) 3 / 2 

and Q is the stellar rotation rate seen by a distant inertial observer. 

At periastron, we have Q < 49.4 - 113.9 = -64.5 rad/yr. This can be compared with 
the critical rotation rate of star A: fi c = a/ GMa/ R\ = 60.0 rad/yr. This means that, in 
absolute terms, the star would have to be rotating faster than critical. This problem can 
be reduced by choosing n > 2. But the larger n is, the more d9cm/dt will be dominated 
by orbital rotation, which is highly non-uniform around the orbit, removing the need for 
variable detectability (see below). There are other problems with this model, however. For 
example for small n, Q is in a direction opposed to that of the orbit. In addition this model 
only works if the orbital and rotational plane of star A are co-planer, which may be unlikely 
for such a wide binary. 

The rotation of star A will necessarily impose curvature on any CIR originating near 
the surface of the star (e.g. Cranmer & Owocki 1996). However, to explain the observed 
strongly varying time-separation around the orbit requires dv/dt to dominate over \Q\ at 
most phases. Because of this, and for the sake of simplicity, we ignore curvature of the CIRs 
in what follows. Then the variation of the flare FWHM and the average flare separation as 
a function of time can be expressed as follows: 

At(FWHM) = (26 x)l {dv/dt - O) = — ^ — , ' ^rj^—, (4) 

v ' y " y 1 ' tt(1 + ecosi;) 2 - (1 -e 2 ) 3 / 2 Pfi/2 w 

and 

At(sep) = —t^t, ^ = -^-At(FWHM). (5) 

v F> n (dv/dt -n) n6 x y ' y ' 

The CIR model directly predicts a linear correlation between the FWHMs and flare separa- 
tions. Fitting eq. [4] to the observed FWHM vs phase yields 9x = 0.08 radians or 9x = 4° 
with a value of fl = 0.014 days -1 or P rot = 450 days. We then constrained Ox between 0.06 
and 0.1 radians, and fit eq. [5] to the observed FWHM vs. average flare separation. This 
yielded n = 75 co-rotating interaction regions around the star. This value of n is far higher 
than is normally encountered in OB stars, while Ox is relatively small. In a sense, such 
a large value of n more closely matches the notion of a large number of clumps. Further- 
more, allowing for curvature, which is stronger for larger Q, also introduces an asymmetry 
in At(FWHM) or At(sep) vs. phase, depending on whether one is receding from or ap- 
proaching periastron, which is not seen. Overall, then, it appears difficult to explain 7] Car's 
X-ray flares via CIR action. 
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4.2. Flares from Intrinsic Instabilities in the Colliding-wind Region 

Numerical smoothed particle hydrodynamics (SPH) models of the X-ray emission pro- 



duce by the wind- wind collision in rj Carinae by Okazaki et al. (2008) shows "flare-like" 
instabilities especially prominent just before periastron. The behavior of the SPH mod- 
els is qualitatively similar to the observed behavior of the RXTE X-ray lightcurve. These 
variations in the model lightcurves were due to transient absorption of the X-ray emission 
(assumed to originate from a point-like region near the stagnation point) by blobs produced 
near the unstable wind- wind interaction surface passing through the observer's line of sight. 
As discussed above the X-ray spectra from a number of flares shows that the flares instead re- 
sult from excess high energy emission rather than increases in absorption. But the wind-wind 



interaction in r\ Carinae may be radiatively unstable near periastron (Pittard & Corcoran 



2002 Parkin et al.||2009 ), and it's worth considering whether such instabilities could produce 



the flare-like behavior we observe. 

Any gas with density p has a preferred timescale for oscillations, which goes as the free- 
fall time: St ~ tff = 0.54/(Gp) 1 / 2 . Assuming that the density in the bow-shock head where 
the excess X-rays arise is proportional to the arriving wind density at the point of impact, 
and with M = Anr 2 p{r)V (r) , one finds St oc rV^r) 1 ^ 2 . Thus, for V(r) ~ and constant 
for most of the orbit, St oc r, as observed for At(FWHM). All this may explain the time 
scale variations (widths of flares). But density variations produced by stochastic instabilities 
should show both overdensities and underdensities, which would produced bright peaks and 
fainter troughs, respectively. While we see clear evidence of brightness enhancements, there 
is little evidence for such troughs in the RXTE X-ray lightcurve, especially away from the 
X-ray minimum in the interval 0.2 < < 0.8, where there is less crowding and troughs 
should be more obvious if they existed. Also, the observed excess peaks are much narrower 
than any possible deficit, contrary to what is expected if due to oscillations. Because of 
these difficulties, the whole idea of oscillations does not seem too convincing, although more 
detailed considerations may be useful. Such considerations require detailed numerical models 
which are beyond the scope of this paper. 



5. Conclusions 

Basically, all three models have varying degrees of problems. In addition, there is the 
shared problem that the orbital parameters of r) Carinae are not very well known. But despite 
the inherent subjective problem of actually defining the flares in t] Carinae we believe that 
the clump model is the simplest (and therefore most likely) explanation yet found for the 
flares. LBV stellar pulsations on a time scale of ~85 d seem less likely because of the lack 
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of a true oscillation signature. 



Clumping on relatively large scales has been invoked to explain the X-ray lines in O 



stars (Feldmeier et al. 2003). Walter & Zurita Heras (2007) use similar model parmeters 



to explain flaring in X-ray binaries. Thus, wind collisions have some similarity to massive 
X-ray binaries, in which relatively hard X-ray flares are seen as a likely result of clumpy 
winds being accreted by the compact component. If so, then one should should see X-ray 
flares in other wind-wind colliding systems, such as the prototype WR + O system WR140 
= HD 193793, WC7pd + 05, P = 7.94 yr, e = 0.88. The X-ray lightcurve of WR140 is also 
being monitored by RXTE. If it shows flares (in fact there is some evidence for absorption 
dips possibly from clumps), then they are far less obvious than those in 77 Carinae. One 
difference between these two systems, though, is that the two colliding winds in WR140 are 
both very fast and may not contain the prominent clumps that LBV winds appear to have 



(Davies et al. 2007), leading to the giant X-ray flares that make 77 Carinae so spectacular. 
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Table 2. Flare Parameters 
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Table 2 — Continued 
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a Some gaps or additions occur in the flare numbering due to omission after initial as- 
signment too close to the noise level or two originally selected flares (nos. 7 & 30) being 
subsequently split into two for greater consistency, while no. 8 was left as one. 

b Flares are not observed during the X-ray minima, so we do not calculate the average 
separation for the last flare before the minimum or the first flare after minimum. 
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Fig. 1. — Top: RXTE X-ray lightcurve (solid black line), adopted brightness level for un- 
derlying (baseline) X-ray variations (dashed line), and residuals (solid gray line). The filled 
circles show the points adopted to represent the underlying brightness due to the colliding- 
wind emission. Tick marks indicate the positions of the peaks of the identified flares, along 
with the flare identifications from Table [2] (except just prior to minima when the flare spac- 
ing is too small). Bottom: Zoomed plots of the strong flare intervals around the 3 X-ray 
minima; tick marks show flare peaks along with the flare ID numbers from Table [2] 
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Fig. 2. — Top: Comparison of residual net rates (after subtracting the adopted "baseline") 
for the three cycles seen by RXTE. Green is cycle 1, blue cycle 2, and red cycle 3. Bottom: 
Close-up of phase interval 0.9 < < 1.0. 
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Fig. 3. — Full width half maximum (in days) for the identified flares listed in Table |2] vs. 
orbital phase. Green symbols are from cycle 1, blue symbols cycle 2, and red symbols cycle 
3. The smooth curves are the best-fit models described in the text: clump model (see eq. 
[2]), long-dashed line; CIR model, short dashed line. 
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Fig. 4. — Flare FWHM vs. average separation between successive flares. Green symbols are 
from cycle 1, blue symbols cycle 2, and red symbols cycle 3, along with the best clump-model 
fit and the best fit CIR model. 
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Fig. 5. — Plot of flare strength (peak x FWHM) vs. phase. 
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Fig. 6. — Plot in residual intensity units (counts/s) of cycle 1 (green squares) and cycle 3 
(red circles) vs. cycle 2. We interpolated the phasing of the cycle 3 curve to that of cycle 2 
before plotting. In order to avoid spurious artifacts, we interpolated the cycle 2 data to the 
(poorer-sampled) cycle 1 phasing before plotting. The solid circle is centered on (0,0) with 
a radius of 2 counts/s, below which random errors dominate. 
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Fig. 7. — Histogram of residual flare intensity ( i.e. all intensity points after subtracting 
baseline). Turnover occurs at an intensity of about 1.3 cts s -1 . Above this value the dis- 
tribution is characterized by a power-law of slope —2.86 ± 0.31 as shown by the dashed 
line. 
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Fig. 8. — Histogram of the distribution in polar angle of the residual intensity (see figure 
[61). The horizontal line is a constant mean value between 0° and 90°, which is an adequate 
description of the distribution in this angular range (xl = 1-01). 
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Fig. 9. — Discrete Fourier Transform of residuals (amplitude vs. frequency). Low-frequency 
amplitude is shown in the inset. 
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Fig. 10 — Binned CHANDRA MEG +1 order X-ray spectrum of the peak of Flare #43 
(black) from Dec 8, 2008 compared to the X-ray spectrum 4 days after the peak, near flare 
minimum (red). The solid lines are absorbed APEC thermal spectra fit to each spectrum 
individually. Note that the single temperature fit does not account for the helium-like Si XIII 
triplet near 1.85 keV or the S XV triplet near 2.45 keV which are prominent in the spectra. 
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Fig. 11. — Sketch showing a cross section in the orbital plane with idealized, outward propa- 
gating, homologously expanding clumps of increased radial extent Ar c i and of angular width 
A8 m 12° as in Davies et al. (2007) at three different distances along the same trajectory 
from star A (left), along with the companion (star B) and the associated shock cone, at 
two different separations. The narrow region of angular width 28x as seen by star A, where 
wind-collision X-rays arise, is indicated by a zig-zag line at the bow-shock head. Note that 



9x « 0, since X-rays are produced mostly where the winds collide nearly head on (Henley 



et al. 2008). Not to scale. 
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Fig. 12. — Cartoon showing a cross section in the orbital plane with four idealized CIRs 
rotating counterclockwise with rj Carinae A (left). One of the CIRs is shown entering (or 
being entered by - depending on the relative angular speed of each) the narrow region at the 
orbiting bow shock head (angular size 29x) where X-rays arise. Not to scale. 



